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Abstract 

We  evaluate  the  k  —  e  turbulence  model  using  di¬ 
rect  numerical  simulation  (DNS)  data  of  a  Mach  4 
boundary  layer.  We  find  that  the  low  Reynolds  num¬ 
ber  damping  functions  for  the  Reynolds  stress  must 
be  corrected  by  the  density  ratio  to  match  the  DNS 
data.  We  present  the  budget  of  the  k  equation  and 
assess  the  modeling  of  the  various  source  terms.  The 
models  for  all  the  source  terms,  except  for  the  produc¬ 
tion  and  dilatational  dissipation  terms,  are  found  to 
be  adequate.  Finally,  we  present  the  solenoidal  dis¬ 
sipation  rate  equation  and  compute  its  budget  using 
the  boundary  layer  data.  We  compare  this  equation 
with  the  dissipation  rate  equation  in  an  incompressible 
flow  to  show  the  equivalence  between  the  two  equa¬ 
tions.  This  is  the  basis  for  modeling  the  solenoidal 
dissipation  equation.  However,  an  additional  term  in 
the  equation  due  to  variation  of  fluid  viscosity  needs 
to  be  modeled. 

1.  Introduction 

The  k  —  e  turbulence  model  is  widely  used  to  sim¬ 
ulate  incompressible  and  compressible  flows.  In  this 
model,  transport  equations  for  the  turbulent  kinetic 
energy,  k,  and  its  dissipation  rate,  e,  are  solved. 
The  Reynolds  stress  is  then  modeled  in  terms  of  k 
and  e,  along  with  a  damping  function  to  account 
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for  the  low  Reynolds  number,  i?e,  effects  close  to  a 
solid  wall.  The  modeling  of  the  unclosed  terms  in  k 
and  e  equations,  and  the  low  Re  damping  function  are 
mostly  based  on  dimensional  arguments.  The  validity 
of  these  assumptions  often  limit  the  performance  of 
the  model  when  applied  to  engineering  problems.  The 
k  —  e  turbulence  model  has  been  tested  against  a  wide 
range  of  experimental  data.  However,  most  of  the  data 
are  limited  to  the  mean  flow  quantities.  The  higher 
order  correlations  involved  in  the  unclosed  terms  are 
difficult  to  measure  experimentally.  By  comparing  the 
model  prediction  of  the  mean  flow  quantities  with  the 
experimental  data,  one  can  assess  the  overall  perfor¬ 
mance  of  the  turbulence  model  but  cannot  evaluate 
the  assumptions  made  for  each  unclosed  term.  In  this 
regard,  a  direct  numerical  simulation  (DNS)  database 
is  very  useful,  wherein  the  unclosed  terms  can  be  eval¬ 
uated  exactly  and  compared  to  their  modeled  counter¬ 
part. 

DNS  of  simple  flows  have  been  previously  used  to 
develop  and  test  turbulence  models.  However,  most  of 
the  work  has  been  done  for  incompressible  flows.  Man- 
sour  et  al . 1  used  DNS  data  of  channel  flows  to  com¬ 
pute  the  budget  of  the  Reynolds  stresses  and  the  dis¬ 
sipation  rate  of  the  turbulent  kinetic  energy.  The  ex¬ 
isting  closure  models  for  the  source  terms  were  tested 
against  the  DNS  data  and  were  found  to  be  inade¬ 
quate.  Rodi  and  Mansour2  evaluated  various  low  Re 
damping  functions  against  DNS  data  of  channel  and 
boundary  layer  flows.  They  also  tested  existing  models 
for  the  e-equation  and  proposed  new  models  based  on 
the  DNS  data.  Also,  Nagano  and  Shimada3  presented 
comprehensive  DNS-based  modeling  for  the  Reynolds 
stress,  and  the  source  terms  in  the  k  and  e  equations. 
The  latest  work  in  this  direction  is  by  Sarkar  and  So4 
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Model 

Code 

D 

fu. 

Launder-  Sharma 

LS 

2  V\fk:y 

exp  (—3.4/ (1  +  0.002ReT)2) 

Chien 

CH 

2 vk/y2 

1  —  exp(— 0.0115y+) 

Lam-Bremhorst 

LB 

0 

[1  -  exp  (— 0.0165y*)]2  (1  +  20.5 /ReT) 

Shih-Mansour 

SM 

2V'/k:y2 

1  -  exp(— 6  x  10-3y+  -  4  x  10~4y+2 

+2.5  x  10~  V"3  -  4  x  10~9y+4) 

N  agano-Tagawa 

NT 

0 

[1  -  exp  (-y+/26)]2  (l  +  4.1/Re^4) 

Rodi-Mansour 

RM 

2  V^fk^y 

1  -  exp(— 2  x  10 -Ay+  -  6.5  x  10 ~V“2) 

Myong-Kasagi 

MK 

0 

(1  +  3.45 /y/Rer)  [1  -  exp(-y+/70)J 

Table  1.  Damping  function,  /M,  used  in  different  low  Re  versions  of  the  k  —  e  turbulence  model.  Here, 
y+  =uTy/vw,  y*  =  k1/2y!vw  and  ReT  =  k2/ve. 


where  DNS  data  for  complex  flows  is  used  for  evaluat¬ 
ing  models. 

DNS  based  modeling  for  compressible  flows  is  more 
limited  than  for  incompressible  flows.  The  transport 
equation  for  the  turbulent  kinetic  energy  in  a  com¬ 
pressible  flow  is  similar  to  the  incompressible  equation, 
but  has  additional  source  terms  due  to  compressibility. 
Huang  et  al .6  use  channel  flow  data  at  Mach  1.5  and 
3.0  to  compute  the  budget  of  these  source  terms.  They 
present  detailed  analysis  of  the  compressibility  terms 
and  develop  models  for  them.  Similar  results  are  pre¬ 
sented  by  Guar  ini  et  al.6  for  a  Mach  2.5  boundary 
layer.  However,  in  contrast  to  incompressible  flows, 
there  is  no  comprehensive  testing  of  the  models  for  all 
the  source  terms  in  the  fc-equation. 

The  usual  approach  for  computing  the  dissipation 
rate  in  a  compressible  flow  is  to  solve  a  modeled  trans¬ 
port  equation  for  the  solenoidal  part,  es,  and  model 
the  dilatational  part,  erl ,  in  terms  of  the  former.  It 
is  assumed  that  es  follows  the  same  dynamics  as  e 
in  an  incompressible  flow,  and  therefore  the  modeled 
e-equation  for  incompressible  flows  is  used  for  comput¬ 
ing  es .  This  approach  is  universally  accepted,  however 
there  is  very  little  understanding  of  the  underlying  as¬ 
sumption.  Specifically,  there  is  no  rigorous  study  of 
how  the  exact  form  of  the  transport  equation  for  es 
relates  to  the  incompressible  e-equation,  k  and  e  are 
computed  from  the  modeled  transport  equations  and 
then  used  to  compute  the  Reynolds  stresses.  The  low 
Re  damping  functions  developed  for  incompressible 
flows  are  used  in  the  compressible  case.  Once  again, 
there  has  been  no  evaluation  of  these  damping  func¬ 
tions  against  compressible  DNS  data.  To  summarize, 


more  research  is  required  to  test  the  k  —  e  model  for 
its  use  in  compressible  flows.  Evaluating  the  model 
against  DNS  data  will  lead  to  better  understanding  of 
its  limitations  so  that  it  can  be  improved  to  predict 
compressible  turbulent  flows  more  reliably. 

In  this  paper,  we  assess  the  accuracy  of  the  k  —  e 
turbulence  model  using  the  DNS  database  of  a  Mach 
4,  adiabatic,  turbulent  boundary  layer.7  In  Section  2, 
we  test  different  low  Re  versions  of  the  k  —  e  model 
for  the  Reynolds  stress.  In  Section  3,  we  present  the 
budget  of  the  turbulent  kinetic  energy  transport  equa¬ 
tion  and  evaluate  existing  models  for  the  source  terms 
against  the  DNS  data.  In  Section  4,  we  present  a  simi¬ 
lar  budget  of  the  transport  equation  for  the  solenoidal 
dissipation  and  address  its  modeling.  We  study  how 
this  equation  relates  to  the  incompressible  dissipation 
equation.  As  pointed  out  earlier,  this  is  the  basis  for 
modeling  the  es -equation.  The  similarities  and  the 
differences  between  the  two  equations  are  pointed  out, 
followed  by  their  modeling  implications.  Conclusions 
are  presented  in  Section  5. 

2.  Reynolds  Stress 

The  Reynolds  stress  tensor,  ,  is  defined  as, 

pTij  =  pu’’u’’  =  p  u”u”  (1) 

where  the  overbar  and  tilde  represent  Reynolds  and 
Favre  averaging,  respectively,  u"  is  the  Favre-fluctuating 
velocity  component  and  p  is  the  density.  In  a  zero 
pressure  gradient  boundary  layer,  only  Txy  =  u"v"  is 
important,  where  u  and  v  are  the  velocity  compo¬ 
nents  in  the  streamwise,  x,  and  wall  normal,  y,  direc- 
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Figure  1.  Normalized  Reynolds  stress  in  a  Mach  4 
boundary  layer  obtained  using  DNS. 


Figure  2.  Variation  of 
boundary  layer. 


fjtNS  with  V+  in  a  Mach  4 


tions.  Figure  1  shows  the  variation  of  rxy  normalized 
by  u2  across  the  boundary  layer  as  computed  using 
the  DNS  data.  uT  is  the  friction  velocity  defined  as 
pwv%.  —  tw  ,  where  pw  and  tw  are  the  density  and  the 
shear  stress  at  the  wall,  respectively.  Guarini  et  al.6 
present  similar  data  for  a  Mach  2.5  boundary  layer. 
They  have  shown  that  the  normalized  Reynolds  stress 
in  the  compressible  flow  matches  the  incompressible 
data  by  Spalart8  when  scaled  by  the  density  ratio, 
p/ Pw  ■  We  get  similar  results  using  the  Mach  4  bound¬ 
ary  layer  data. 

As  per  the  Boussinesq  eddy-viscosity  approxima¬ 
tion,  u"v"  in  a  zero  pressure  gradient  boundary  layer 
is  modeled  as 

-pu"v"=pT^-  (2) 


The  eddy  viscosity,  pt  ,  is  modeled  in  terms  of  k  and 


,  Pk2 

Pt  —  Jf, - , 


(3) 


where  cy  =  0.09  is  a  model  constant  and  /;t  is  a  damp¬ 
ing  function.  /p  =  1  away  from  the  wall  and  it  is  less 
than  1  close  to  the  wall  to  account  for  the  effect  of 
viscosity  and  the  damping  of  turbulence  by  the  solid 
wall.  Different  low  Re  versions  of  the  k  —  e  turbulence 
model  use  different  forms  of  /M.  Rodi  and  Mansour2 
evaluated  some  commonly  used  models  against  DNS 
data  for  incompressible  channel  and  boundary  layer 
flows.  They  compare  the  model  damping  functions  to 
an  equivalent  computed  from  the  DNS  data  as  fol¬ 
lows, 


J  fl 


cyk2  du/dy 


(4) 


where  Reynolds  averaging  inRef.  2  is  replaced  by  Favre 
averaging  for  compressible  flows.  Fig.  2  shows  the  vari¬ 
ation  of  /°NS  with  y+  as  computed  using  DNS  data 
of  the  Mach  4  boundary  layer.  Here,  y+  =  yuT/V w  is 
the  non-dimensional  distance  from  the  wall,  u  is  the 
kinematic  viscosity  and  the  subscript  w  represents  the 
value  at  the  wall.  As  pointed  out  earlier,  =  1  away 
from  the  wall  corresponds  to  =  0.09.  However,  in 
the  boundary  layer  data  in  Ref.  2  ,  /°NS  attains  a 
value  lower  than  1  away  from  the  wall  (/^NS  =  0.83 
at  y+  =  100).  Thus,  the  value  of  is  0.075  so  that 
^dns  ^  awa y  from  the  wall.  In  the  Mach  4  bound¬ 
ary  layer  data  used  here,  /°NS  attains  a  value  of  0.40 
at  about  y+  =  100.  This  is  much  lower  than  the  value 
reported  inRef.  2 for  the  incompressible  case.  However, 
similar  to  the  Reynolds  stress,  scaling  f fl  with  the  den¬ 
sity  ratio  improves  the  comparison  with  the  incom¬ 
pressible  damping  function.  /#l  =  (p/pw)fn  is  higher 
than  (Fig.  2)  and  is  about  0.70  at  y+  =  100,  which 
is  much  closer  to  the  incompressible  value.  Also,  using 

=  0.07,  instead  of  the  traditional  value  of  0.09,  in 
the  compressible  boundary  layer  we  get  ~  1  away 
from  the  wall. 

Next,  we  compare  different  low  Re  k  —  e  models  for 
the  Reynolds  stress  against  the  boundary  layer  data. 
The  damping  functions  for  the  models  by  Launder  and 
Sharma  (LS),9  Chien  (CH),10  Lam  and  Bremhorst 
(LB),11  Shih  and  Mansour  (SM),12  Nagano  and 
Tagawa  (NT),13  Rodi  and  Mansour  (RM),2  and  My- 
ong  and  Kasagi  (MK)14  are  given  in  Table  1.  It  is  to 
be  noted  that  some  models  use  a  modified  dissipation 
rate,  e,  in  place  of  e  in  Eq.  (3).  e  =  e  —  D ,  where  D 
for  each  model  is  given  in  Table  1.  To  facilitate  com- 
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Figure  3.  Comparison  of  model  damping  functions 
against  /M  computed  from  DNS  data  of  a  Mach  4 
boundary  layer.  The  model  codes  are  given  in  Table  1. 


Figure  4.  Budget  of  fe-equation  using  DNS  data  of 
a  Mach  4  boundary  layer.  The  source  terms  are  nor¬ 
malized  by  pwu^/vw. 


parison  with  the  DNS  data,  we  compute  a  modified 
form  of  the  model  damping  functions,  =  /Me/e  and 
evaluate  it  against  /^NS  (obtained  using  =  0.07). 
The  results  are  presented  in  Fig.  3.  We  see  that  all 
the  models  over-predict  the  DNS  data  in  the  region 
y+  =  10  through  40,  but  the  final  approach  of  /J)NS  to 
1  is  better  predicted  by  the  model  by  SM.  On  the  other 
hand,  the  models  NT  and  RM  approach  1  faster  than 
/^NS ,  whereas  the  models  LB  and  MK  approach  1  rel¬ 
atively  slow.  Some  of  these  models  are  tested  against 
data  for  an  incompressible  boundary  layer  by  Rodi  and 
Mansour2  and  the  comparison  shown  here  is  similar 
to  their  results.  The  behavior  of  the  models  LS  and 
CH,  not  shown  here,  is  also  similar  to  the  incompress¬ 
ible  case.  The  damping  function  for  LS  rises  to  1  far 
too  quickly  whereas  f fl  for  CH  is  much  lower  than  the 
DNS  data  and  approaches  1  very  slowly.  Overall,  the 
model  by  SM  is  found  to  be  closest  to  the  DNS  data. 


3.  Turbulent  Kinetic  Energy 

In  this  section,  we  present  a  budget  of  the  turbulent 
kinetic  energy  computed  using  the  boundary  layer  data 
and  evaluate  models  for  each  source  term. 


Budget  of  transvort  equation 
The  turbulent  kinetic  energy  in  a  compressible  flow 
is  defined  as 

fc=|«,  (5) 

and  its  transport  equation  is  given  by, 6 
Dpk 


Dt 


—  Pk  +  Tfc  +  Dk  —  pe  +  nfc  +  C\  +  C2  +  O3  .  (6) 


The  left-hand  side  represents  the  substantial  deriva¬ 
tive  of  pk  and  it  balances  the  source  terms.  The  first 
five  source  terms  are  turbulent  production,  turbulent 
transport,  viscous  diffusion,  viscous  dissipation  and 
pressure  transport  of  turbulent  kinetic  energy,  respec¬ 
tively.  The  last  three  terms  are  due  to  compressibility. 
The  expressions  for  these  source  terms  are  given  below, 


Pk 

Tk 

Dk 


,dui 


-  (k!«) 


<3U'id 

Hfc  =  -  (pX)  . 
Cl  =  -Cjp.i 

Cl  —  ui  &ij,j 

C3  =  p'Uii 


where,  p  is  the  pressure,  Oij  is  the  viscous  stress  ten¬ 
sor  and  the  prime  represents  the  Reynolds-fluctuating 
quantities. 

In  a  zero  pressure  gradient  boundary  layer,  the  un¬ 
steady  and  convection  terms  on  the  left-hand  side  of 
Eq.  (6)  are  negligible  and  the  turbulent  kinetic  en¬ 
ergy  is  mainly  governed  by  the  source  terms.  We  com¬ 
pute  the  magnitude  of  these  source  terms  using  the 
boundary  layer  data  and  show  the  budget  in  Fig.  4. 
The  source  terms  are  normalized  by  pwu^/Vw  .  We  see 
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Figure  5.  Comparison  of  models  for  the  production 
of  turbulent  kinetic  energy,  Pk ,  against  DNS  data  of 
a  Mach  4  boundary  layer. 


Figure  6.  Comparison  of  models  for  the  turbulent 
transport,  Tk ,  and  pressure  transport,  11^ ,  terms 
against  DNS  data  of  a  Mach  4  boundary  layer. 


that  the  budget  is  very  similar  to  the  data  presented 
by  Guarini  et  al.6  for  a  Mach  2.5  boundary  layer.  In 
the  high  Reynolds  number  region,  y+  >  50,  produc¬ 
tion,  Pk ,  and  dissipation,  pt ,  are  dominant  and  the 
other  terms  are  comparatively  small.  Closer  to  the 
wall,  production  and  dissipation  are  still  the  largest 
terms  but  turbulent  transport,  Tk ,  and  viscous  diffu¬ 
sion,  Dk ,  are  also  important.  The  viscous  diffusion 
reaches  a  maximum  at  the  wall  and  balances  the  non¬ 
zero  dissipation  rate.  The  rest  of  the  terms  are  zero  at 
the  wall.  The  pressure  diffusion  term  attains  a  max¬ 
imum  value  of  0.025  near  the  wall  and  is  otherwise 
small  throughout  the  boundary  layer.  The  sum  of  the 
three  compressibility  terms,  C ,  is  also  shown.  This 
term  has  been  found  to  be  negligible  in  a  Mach  2.5 
boundary  layer  by  Guarini  et  al.6  and  in  Mach  1.5 
and  Mach  3.0  channel  flows  by  Huang  et  al.5  In  our 
boundary  layer  data  at  Mach  4,  we  find  that  the  max¬ 
imum  value  of  C  is  0.013  at  about  y+  =  12 .  Also,  the 
magnitude  of  the  compressibility  terms  C\  and  C3  are 
negligible  compared  to  Ci . 


Modeling  of  source  terms 

In  a  boundary  layer,  the  production  term  simplifies 
to 

—  3u 

Pk~-pu"v"—.  (7) 

ay 

The  modeling  of  the  Reynolds  stress,  u"v" ,  in  terms 
of  the  turbulent  eddy  viscosity  is  discussed  earlier.  We 
use  some  of  the  models  for  eddy  viscosity  presented  in 
section  2  to  compute  Pk  and  the  results  are  shown  in 


Fig.  5  along  with  the  value  of  Pk  as  obtained  from 
the  DNS  data.  We  see  that  the  models  RM  and  NT 
reproduce  the  initial  increase  in  Pk  but  over-predict 
the  peak  value  by  factors  of  2  and  3,  respectively.  The 
prediction  by  SM  is  in  between  RM  and  NT.  The  re¬ 
maining  models,  not  shown  in  Fig.  5,  predict  even 
higher  values  of  the  peak  production.  It  is  to  be  noted 
that  all  the  models  predict  a  maximum  in  Pk  at  about 
y+  =  10  which  matches  the  location  of  the  peak  in  the 
Pk  profile  computed  from  the  boundary  layer  data. 
Also  the  models  predict  Pk  in  the  region  y+  >  40 
fairly  well. 

Next,  we  consider  the  turbulent  transport,  Tk ,  and 
pressure  transport,  ,  terms.  These  two  terms  are 
modeled  together  in  terms  of  the  turbulent  viscosity 
and  the  gradient  of  the  turbulent  kinetic  energy, 


Tk  +  nfc 


f  f^rdk\ 

\ak  dyj 


(8) 


Here,  ak  =  1.3  is  the  Prandtl  number  correspond¬ 
ing  to  the  turbulent  diffusion  of  k.  To  check  the  va¬ 
lidity  of  the  gradient  transport  assumption,  we  first 
compute  ht  using  DNS  data  for  u"v"  and  du/dy  in 
Eq.  (2).  Using  this  value  of  the  eddy  viscosity,  denoted 
by  /i°NS ,  we  can  eliminate  the  errors  due  to  modeling 
of  ht  ■  We  see  that  {tixNSky)y/crk  reproduces  the  vari¬ 
ation  in  Tk  +  nfc  with  y  but  under-predicts  the  mag¬ 
nitude  of  the  peaks  (Fig.  6).  Mansour  et  al.15  show 
similar  comparison  between  the  model  and  the  source 
term  in  an  incompressible  channel  flow.  This  may  be 
due  to  too  high  a  value  of  ak .  Next,  we  compute 
Ht  using  the  turbulence  models  discussed  earlier.  We 


Figure  7.  Comparison  of  the  gradient  transport 
model  for  the  viscous  diffusion  term,  Dk ,  against  DNS 
data  of  a  Mach  4  boundary  layer. 


have  seen  in  Fig.  2  that  the  models  over-predict  the 
value  of  in  the  range  y+  <  25.  These  modeled 
values  of  pt  result  in  a  higher  value  of  the  peaks  of 
the  model  expression  in  Eq.  (8),  which  match  Tk  +11*, 
better  than  /i^NS .  The  model  by  RM  is  closest  to 
the  DNS  data,  whereas  other  models  over-predict  the 
magnitude  of  T*  +  IR .  The  viscous  diffusion  term, 
Dk ,  is  modeled  in  a  manner  similar  to  the  turbulent 
transport  term  using  the  gradient  transport  hypothe¬ 
sis.  Huang  et  al 5  have  presented  a  detailed  analysis 
of  this  term  to  show  that 


Dk 


(9) 


Figure  8.  Total  dissipation  rate,  e,  and  its  solenoidal, 
es ,  and  dilatational,  ed,  parts  as  computed  using  DNS 
data  of  a  Mach  4  boundary  layer.  All  the  dissipation 
rates  are  normalized  by  uf  / uw .  aiMf  is  also  pre¬ 
sented  for  comparison. 


et  al.6  The  ratio,  ed/es  is  less  than  3%  throughout 
the  boundary  layer.  Huang  et  al.  have  also  reported 
that  es/e  is  as  low  as  0.8  in  the  viscous  sub-layer.  This 
is  largely  due  to  the  fact  that  correlation  between  the 
viscosity  fluctuations  and  the  velocity  gradient  con¬ 
tributes  to  the  total  dissipation.  However,  in  our  adi¬ 
abatic  boundary  layer,  we  find  that  es/e  ~  1  (Fig.  8). 

The  solenoidal  dissipation-rate  is  computed  from  a 
modeled  transport  equation  which  is  discussed  in  de¬ 
tail  in  the  next  section.  The  dilatational  dissipation  is 
modeled  in  terms  of  es  as 


where  p  is  the  mean  molecular  viscosity  of  the  fluid. 
We  see  that  the  above  equation  is  approximately  true 
for  this  test  case  (Fig.  7). 

The  dissipation  rate,  e,  of  the  turbulent  kinetic  en¬ 
ergy  can  be  split  into  the  solenoidal  part,  es ,  the  di¬ 
latational  part,  ed,  the  anisotropy  part  and  additional 
terms  due  to  fluctuations  in  the  molecular  viscosity. 5 


pes  =  p  lj'm'a  with  w  =  V  x  u , 

_L _  (io) 

ped  =  |  p  u'^u'k  k  , 

where,  w  is  the  vorticity.  Fig.  8  shows  the  variation 
of  the  total  dissipation,  the  solenoidal  dissipation  and 
the  dilatational  dissipation  in  the  boundary  layer.  All 
the  quantities  are  normalized  by  u\/vw.  We  see  that 
ed  “C  es  which  is  similar  to  the  results  obtained  by 
Huang  et  al.5  for  the  channel  flows  with  cooled  walls 
and  for  a  boundary  layer  on  adiabatic  wall  by  Guarini 


ed=onM?es,  (11) 

where  Mt  =  \/2 k/a  is  the  turbulent  Mach  number,  a 
is  the  Favre-averaged  speed  of  sound  and  cti  =1.0  is 
the  generally  accepted  model  constant.  We  see  that 
ot\ Mf ,  as  computed  using  the  DNS  data,  is  much 
larger  than  ed/es .  Thus,  the  above  model  greatly  over¬ 
predicts  the  dilatational  dissipation.  This  is  similar  to 
the  observations  in  other  compressible  flows.5,6 

The  compressible  source  terms  C\  and  C2  involve 
u "  which  needs  to  be  modeled.  The  values  of  the 
stream- wise,  u" ,  the  wall  normal,  v" ,  and  the  span- 
wise,  w” .  fluctuations  normalized  by  the  edge  veloc¬ 
ity,  Ue ,  are  presented  in  Fig.  9.  We  see  that  u"  is  the 
largest  of  the  three  components  and  attains  a  peak 
value  of  about  15%  of  the  edge  velocity.  The  wall- 
normal  component  is  about  2%  of  Ue,  whereas  the 
span-wise  component  is  almost  negligible  throughout 
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Figure  9.  Comparison  of  models  for  u" ,  v"  and  w" 
against  DNS  data  of  a  Mach  4  boundary  layer.  The 
two  models  are  HCB  by  Huang  et  al.5  and  Rubesin.16 


the  boundary  layer.  Huang  et  al.5  proposed  the  fol¬ 
lowing  model  for  u" , 


1  u"T" 


u ■  = 


u"v" 


(12) 


where  T  is  the  temperature.  The  model  is  based  on 
compressible  channel  flow  data  where  the  pressure  fluc¬ 
tuations  are  negligible  compared  to  the  density  and 
temperature  fluctuations  (Coleman  et  al.17).  In  the 
Mach  4  boundary  layer  data  used  in  this  paper,  the 
pressure  fluctuations  are  found  to  be  small  and  thus 
the  above  model  is  expected  to  do  well  for  this  test 
case.  Comparing  the  model  against  the  DNS  data  in 
Fig.  9,  we  see  that  the  predictions  for  v"  and  w"  are 
very  good  whereas  u"  is  marginally  over-predicted. 
The  u”  -model  proposed  by  Rubesin16  is  given  by 


k 

e 


dT  1 
~dy  f 


(13) 


where  Ci  is  a  calibration  constant  and  T  is  the  Favre- 
averaged  temperature.  A  value  of  c\  =  0.7  reproduces 
the  peak  in  u"  well  but  over-predicts  the  values  in  rest 
of  the  boundary  layer,  whereas  C2  =0.1  reproduces  v" 
very  well  in  almost  all  of  the  boundary  layer. 


4.  Solenoidal  Dissipation  Rate 


followed  by  the  budget  as  computed  using  the  bound¬ 
ary  layer  data.  The  modeled  es -equation  is  a  direct 
extension  of  the  incompressible  dissipation  equation. 
We  compare  these  two  equations  to  assess  the  accu¬ 
racy  of  this  modeling  approach. 


Budget  of  transport  equation 
The  solenoidal  dissipation  rate,  es ,  is  given  by 

es=p  (u'i  ju'i  j  -  u'iju'j  ^  ,  (14) 

where  9  is  the  mean  kinematic  viscosity.  The  trans¬ 
port  equation  for  es  has  the  following  form 

^  =  P1  +P*  +P'f  +P*  +Tf+Df-Yf  +Tfc+Bf+Ct  . 

(15) 

where  the  expressions  for  the  various  source  terms  are 
listed  below. 


Pe 

=  -^2Kj 

u'j,i )  u'k,j  uhk 

P2 

=  -^2K  j 

-  u),i )  K,k  ™k,j 

P? 

=  ~v2K,j 

~  Uj,i)  U'k  V'ijk 

P? 

=  -^2K  j 

~  U’j,i)  U'i,kU'k,j 

Te  =  -*KK  j-UjJKjlk 
D<  =  ~  <A,i)  tkk 

Ye  =  P22{u'iJ-u'jti),ku,iijk 

Te  =  V  {u'ij  ~  W'J  U'tj  U'kk 

Be  =  9  2  {u'ij-u'j  Jpjpj/p2 

C€  =  9  2  (uij  —  Uj  j)  {y'Ui^kkj  +  tpe)  ■ 


2  K* -«$,<).* 


The  first  four  terms  represent  different  production 
mechanisms.  P1  is  the  mixed  production  term,  P2 
is  the  production  by  mean  velocity  gradient,  Pf  is  the 
gradient  production  term  and  Pf  represents  the  tur¬ 
bulent  production  mechanism.  The  turbulent  trans¬ 
port,  viscous  diffusion  and  viscous  dissipation  of  es 
are  given  by  T€ ,  De  and  Yt ,  respectively,  whereas  the 
last  three  terms  in  Eq.  (15)  are  due  to  compressibility 
effects.  The  compressible  transport  term,  Tec,  repre¬ 
sents  the  effect  of  non-zero  dilatation,  Be  is  due  to 
baroclinic  torques,  and  Ce  includes  the  effect  of  vari¬ 
ation  in  fluid  viscosity.  ipe  is  given  by 


We  - 


Oxh 


+ 


dxk  \_ 


P,k  /  i  \  i  2/b« 

\Ui,k  T  ^k,i)  T  o 

P  3  p  ’ 


k  —  e  modeling  of  compressible  flows  solve  a  mod¬ 
eled  transport  equation  for  the  solenoidal  dissipation. 
Here,  we  start  with  the  exact  form  of  the  es  -equation 


A  budget  of  the  source  terms  in  Eq.  (15)  as  com¬ 
puted  using  the  DNS  data  is  shown  in  Fig.  10.  Follow¬ 
ing  the  work  by  Rodi  and  Mansour2  for  incompressible 
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Figure  10.  Budget  of  the  es -equation  computed  us¬ 
ing  DNS  data  of  a  Mach  4  boundary  layer.  All  the 
terms  are  normalized  by  u^/v^. 

flows,  we  normalize  all  the  source  terms  by  w®  / .  We 
see  that  only  the  turbulent  production  term,  Pf ,  and 
the  viscous  dissipation  term,  Yf ,  are  significant  in  the 
high  Re  region  away  from  the  wall,  and  they  approx¬ 
imately  balance  each  other.  Closer  to  the  wall,  P} 
and  P2  are  as  large  as  P4,  whereas  Pf  is  relatively 
small.  This  pattern  of  the  production  terms  is  very 
similar  to  the  data  presented  by  Rodi  and  Mansour.2 
The  dissipation  term  is  large  in  this  near-wall  region 
to  balance  all  the  production  terms,  whereas  the  vis¬ 
cous  diffusion  term,  De ,  balances  the  large  negative 
dissipation  at  the  wall.  De  is  small  compared  to  the 
production  terms  in  the  rest  of  the  boundary  layer. 
The  magnitude  of  the  turbulent  transport  term  is  of 
the  same  order  as  that  of  the  viscous  diffusion  term, 
except  for  very  close  to  the  wall.  The  compressible 
transport  term,  TJf ,  and  the  baroclinic  term,  Be ,  are 


negligible,  whereas  the  source  term  due  to  variation  in 
fluid  viscosity,  Ce ,  is  relatively  large  in  the  near  wall 
region. 

Modeling  of  transport  equation 
The  modeled  transport  equation  for  es -equation  is 
not  derived  from  its  exact  form  presented  above,  rather 
the  modeled  dissipation  rate  equation  for  incompress¬ 
ible  flows  is  directly  used  to  compute  es  in  a  com¬ 
pressible  flow.  This  is  based  on  the  assumption  that 
es  follows  the  same  dynamics  as  in  an  incompress¬ 
ible  flow.  (Here,  the  subscript  i  is  used  to  represent 
the  incompressible  flow  quantities.)  In  other  words,  it 
is  assumed  that  the  transport  equations  for  es  and  ti 
have  the  same  form.  Here,  we  check  the  validity  of 
this  assumption  by  comparing  the  exact  form  of  the 
two  equations  both  analytically  and  with  the  aid  of 
the  DNS  data. 

The  dissipation  rate  in  an  incompressible  flow  is  de¬ 
fined  as2 

ti  =  v  u'ijuY ,  (16) 

and  the  transport  equation  for  e,  is  given  by 

^-  =  P}+P?+P?+P?+Ti+Ili+Di-Yi  ,  (17) 
where 


Pi 

= 

—  V  2  u'i  jll'k  j  Ui,k 

Pi 

= 

~l'2ui,jUi,k'ak,j 

Pf 

= 

—  V  2  u'i  jlh'f.  Uijk'.  , 

Pi 

= 

-v2u'i,ju'i,ku'k,j  ■ 

Ti 

= 

-"K <A,j\k  > 

n* 

= 

-j/2K ,jp',j\}i/p  ■ 

Di 

= 

"2 

Yi 

= 

^2ui,jkui,jk  ■ 

Once  again,  the  subscript  i  distinguishes  the  source 
terms  in  the  incompressible  case  from  those  in  the 
compressible  equation  (15).  The  nomenclature  of  the 
terms  is  the  same  as  that  in  the  es -equation.  P[  ,  Pf  , 
Pf  and  Pf  represent  the  mixed  production,  the  pro¬ 
duction  due  to  mean  velocity  gradient,  the  gradient 
production  term  and  the  turbulent  production  mech¬ 
anism,  respectively.  The  remaining  terms,  7),  , 

Dj  and  Y) ,  represent  turbulent  transport,  pressure 
transport,  viscous  diffusion  and  viscous  dissipation,  re¬ 
spectively.  Comparing  the  es -equation  with  Eq.  (17), 
we  see  that  the  four  production  terms,  the  turbulent 
transport  term,  the  viscous  diffusion  term  and  the  vis¬ 
cous  dissipation  term  are  present  in  both  the  equa¬ 
tions.  These  are  referred  as  the  common  terms  here¬ 
after.  However,  there  are  additional  terms  —  namely, 
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Figure  11.  Comparison  of  the  terms  in  Eq.  (18)  as 
computed  using  the  Mach  4  boundary  layer  data.  All 
the  quantities  are  normalized  by  (w<w<)ro  . 


Ilj  in  the  incompressible  equation  and  the  compress¬ 
ible  terms  in  the  es -equation.  To  check  whether  es 
follows  the  same  dynamics  as  e, ,  we  compare  the  ex¬ 
pressions  for  es  and  .  We  then  compare  their  trans¬ 
port  equations  -  first,  the  common  terms,  and  then 
the  extra  terms. 

The  solenoidal  dissipation  rate  in  Eq.  (14)  can  be 
written  as  es  =  v  ,  where 


d  _  —  (18) 

..tat  i  at 2  v  ' 


=  <j<j  ~  ~  Q^<e'  +  9'2 


Here,  O'  =  ^  is  the  fluctuating  dilatation.  The  last 

three  terms  represent  the  effect  of  inhomogeneity,  the 
velocity-dilatation  correlation  and  compressibility,  re¬ 
spectively.  Fig.  11  shows  the  magnitude  of  the  terms 
in  the  above  equation  in  the  boundary  layer.  As  noted 
earlier,  9'2  <C  .  We  also  see  that  the  second  and 
the  third  terms  in  (18)  have  negligible  contribution  to 
the  enstrophy.  The  last  three  terms  constitute  ubub  , 
and  therefore 


u'j,Kj «  ■ 


(19) 


Thus,  es  can  be  seen  as  composed  of  two  terms,  the 
first  is  identical  to  and  the  second  part  represents 
the  effect  of  inhomogeneity  and  compressibility.  As 
seen  above,  the  second  term  is  negligible  compared  to 
the  first,  and  therefore 


(20) 


Thus,  es  is  equivalent  to  in  a  compressible  flow. 


Figure  12.  Source  terms  in  the  es -equation  along 
with  the  corresponding  incompressible  terms  (sym¬ 
bols)  given  in  Eq.  (21)  as  computed  using  the  DNS 
data. 


Next,  we  extend  the  above  analysis  for  es  to  the 
common  source  terms  in  the  transport  equation  (15) 
which  have  a  corresponding  term  in  the  -equation. 
Each  source  term  is  split  into  two  parts  -  the  first 
part  is  identical  to  the  respective  incompressible  source 
term  and  the  second  part  represents  the  effect  of  inho¬ 
mogeneity  and  compressibility. 


P!  =  pl  +  2j7uWu'fcj«i,fc 

P?  =  Pi  +29u'j,iu'i,k^j 

p f  =  pf  +  2 


P?  =  Pt  +  2Vu'j^ku'Kj 


p  =  P  +DK  u'iA,j\k 
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We  compute  the  above  terms  using  the  boundary  layer 
data  and  find  that  the  second  part  of  each  source  term 
is  negligible  compared  to  the  first  part.  This  is  analo¬ 
gous  to  the  result  presented  in  Eq.  (19).  Fig.  12  shows 
this  result  where  each  of  the  above  source  terms  is  iden¬ 
tical  to  its  respective  first  part.  Thus,  the  transport 
equation  for  solenoidal  dissipation  simplifies  to 


^  =  Pi  +P?  +P?  +  P*  +  Ti  +  Di  -  Yi+Ct 


In  the  above  equation,  we  have  neglected  the  com¬ 
pressible  transport  term,  Tf ,  and  the  baroclinic  term, 
P€ ,  which  are  shown  to  be  negligible  in  Fig.  10.  Com¬ 
paring  the  simplified  es  -equation  (22)  with  the  incom¬ 
pressible  dissipation  equation  (17),  we  see  that  the  two 
equations  are  identical  except  for  the  terms  11^  and 
Ce .  The  pressure  transport  term,  11^ ,  in  the  incom¬ 
pressible  equation  is  found  to  be  negligible  by  Rodi  and 
Mansour,2  whereas  the  effect  of  the  viscosity  variation 
term,  Ce ,  is  not  negligible  (Fig.  10)  and  should  be  ac¬ 
counted  for.  The  incompressible  models  for  the  rest  of 
the  source  terms  can  be  extended  to  the  compressible 
case. 


5.  Conclusions 

We  evaluate  the  k  —  e  turbulence  model  using  DNS 
data  of  a  Mach  4  boundary  layer.  The  results  are  pre¬ 
sented  in  three  parts.  First,  the  low  Re  models  for  the 
Reynolds  stress  are  tested.  We  find  that  the  models 
must  be  corrected  by  the  density  ratio  to  match  the 
DNS  data.  Next,  we  present  the  budget  of  the  k  equa¬ 
tion  and  test  models  for  each  source  term.  The  mod¬ 
els  predict  too  high  a  peak  for  the  production  term, 
whereas  the  models  for  the  turbulent  transport,  pres¬ 
sure  transport,  viscous  diffusion  and  the  compressible 
terms  perform  reasonably  well.  The  dilatational  dis¬ 
sipation  is  negligible  compared  to  the  solenoidal  part, 
but  the  model  grossly  over-predicts  its  value.  We  com¬ 
pute  the  budget  of  the  solenoidal  dissipation  rate,  es , 
and  find  it  to  be  similar  to  the  budget  of  the  dissipa¬ 
tion  rate,  e« ,  in  an  incompressible  flow.2  The  modeling 
of  the  es  equation  is  based  on  its  equivalence  to  the 
equation.  Comparison  of  the  two  equations,  analyti¬ 
cally  and  with  the  aid  of  DNS  data,  shows  that  they 
are  identical  except  for  a  source  term  due  to  the  vari¬ 
ation  of  fluid  viscosity  in  the  compressible  case.  Thus, 
the  modeled  incompressible  dissipation  rate  equation 
can  be  used  to  compute  es ,  but  the  effect  of  the  addi¬ 
tional  term  must  be  included. 
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